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Abstract 

We calculate a two-pion wave function for the 1 = 2 S-w&ve two-pion system with a finite 
scattering momentum and estimate the interaction range between two pions, which allows us to 
examine the validity of a necessary condition for the finite size formula presented by Rummukainen 

<1^ and Gottlieb. We work in the quenched approximation employing the plaquette gauge action for 

gluons and the improved Wilson action for quarks at 1/a = 1.63 GeV on 32 3 x 120 lattice. The 
quark masses are chosen to give m w = 0.420, 0.488 and 0.587 GeV. We find that the energy 
dependence of the interaction range is small and the necessary condition is satisfied for our range 
of the quark mass and the scattering momentum, k < 0.16 GeV. We also find that the scattering 
phase shift can be obtained with a smaller statistical error from the two-pion wave function than 
qj • from the two-pion time correlator. 
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I. INTRODUCTION 

Calculations of the scattering phase shift represent an important step for expanding 
our understanding of the strong interaction based on lattice QCD to dynamical aspects of 
hadrons. For the simplest case of the 1 = 2 S-w&ve two-pion system, several calculations 
have been reported in Refs. [l|, |2j, y, 14, a, la] • These calculations employ the finite size formula 
presented by Liischer |7| and the extension to the system with a non-zero total momentum 
by Rummukainen and Gottlieb [8j, in which the scattering phase shift is related to the 
energy eigenvalue on a finite volume. In previous applications of the formula the energy was 
calculated from an asymptotic time behavior of the two-pion time correlator. 

The derivation of the finite size formula assumes a condition R < L/2 for the two- 
pion interaction range R and the lattice size L, so that the boundary condition does not 
distort the shape of the two-pion interaction. It is important to examine the validity of 
this necessary condition for reliable results of the scattering phase shift. The CP-PACS 
collaboration calculated the two-pion wave function for the ground state of the 7 = 2 5- 
wave two-pion system and estimated the interaction range from it [9]. In their case the 
scattering momentum is highly small, thus their work is an examination of the condition 
for the scattering length. They found that the interaction range is i? ~ 1.3 — 1.6 fm for the 
pion mass m n = 0.520 — 0.860 GeV. They also evaluated the scattering length and found it 
can be given with a smaller statistical error from the wave function than from the two-pion 
time correlator. 

In the present work we extend their work of the necessary condition for the scattering 
length (the scattering momentum k ~ 0) to that for the scattering phase shift (k ^ 0) 
for the 1 = 2 S-wave two-pion system. In order to set a finite momentum for the ground 
state, we consider a system having the total momentum P = (2n/ L)e x . We calculate the 
two-pion wave function and evaluate the scattering phase shift from it, as carried out by the 
CP-PACS collaboration for the scattering length. 

This article is organized as follows. In Sec. [TT] we briefly review the derivation of the 
finite size formula presented by Rummukainen and Gottlieb |8j, with emphasis on the role 
of the condition for the interaction range. The calculation method of the wave function and 
the simulation parameters are given in Sec. IIII1 In Sec. IIV Al we present our results of the 
wave function and estimate the interaction range. The scattering phase shift from the wave 
function is provided and compared to those with the two-pion time correlator in Sec. HVBl 
Our results of the scattering length and the scattering phase shift at the physical quark mass 
are presented in Sec. IIV CI Our conclusions are given in Sec. |Vj Preliminary reports of the 



present work were presented in Ref. lOj. The calculation was carried out on VPP5000/80 



at the Academic Computing and Communications Center of University of Tsukuba. 



II. FINITE SIZE FORMULA 

We briefly review the derivation of the finite size formula presented by Rummukainen 
and Gottlieb [8], with emphasis on the role of the condition for the interaction range. We 
restrict ourselves the formula for the 1 = 2 S-w&ve two-pion system, but the generalization 
can be easily carried out. The formula has also been derived from another approaches in 






Ref. [ll| and 12|. We follow, however, the original derivation in Ref. 

First we consider a wave function in the infinite volume in the Minkowski space defined 
by 

*(xi',x 2 ) = (0\7r + {x l )n + {x 2 )\7i + {p 1 ),n + {p 2 );m) , (1) 

in the space like region, (xi — x 2 ) 2 < 0. The state \n + (pi), n + (p 2 ); in) is an asymptotic 
two-pion state with the four-dimensional momenta p\ and p 2 . tc + (x) is an interpolating 
operator for ix + at the four-dimensional position x = (x°,x). We assume that the two-pion 
interaction range R is finite and the wave function satisfies 

(D j + m£)tf (xi; X2 ) = o for j = 1, 2 , (2) 

for —(x\ — x 2 ) 2 > R 2 , where Q,- is the d'Alembertian with respect to the coordinate Xj. 

In order to remove a trivial dependence of the center of mass coordinate X = (xi+x 2 )/2, 
we introduce a relative wave function defined by 

<P(x) = y(x i; x 2 )-e iPX , (3) 

where x = X\ — x 2 is the relative coordinate and P = p\ + p 2 = (E, P) is the total four- 
dimensional momentum. From ([2]) <p(x) satisfies 

(□ - k 2 )4>(x) = , (4) 

P ■ d<j)(x) = , (5) 

for —x 2 > R 2 , where k 2 is the scattering momentum defined by 

k 2 = P 2 /A -m\ = (E 2 - P 2 )/4 - ml . (6) 

(Jl} and (JSD also yield 

V 2 - (P • V) 2 /E 2 + k 2 ] <f){x) = . (7) 

We can obtain a relation between the wave function and the scattering phase shift by 
introducing a center of mass frame. The wave function in the center of mass frame 0cm (#cm) 
is related from that in the original frame (f>(x) by the Lorentz transformation, 

0cm(^cm) = (f>(x) . (8) 

xcm is the coordinate in the center of mass frame given by 



CM 



7 (x° - v ■ x) , (9) 

x CM = T[x - vx°] , (10) 



where v is the velocity v = P/E and 7 is the Lorentz boost factor 7 = l/\/I — v 2 = E/yP^. 
The operation 7[x] is defined by 

7 [x] =tx u +x x , (11) 

where x and x ± are components of x parallel and perpendicular to the velocity v, i.e., 
x = v(x ■ v)/v 2 and x ± = x — x . 
From (J4]) and ([5]), 0cm (x) satisfies 

(n-k 2 )<p CM (x) = o, (12) 

^cm(*) = , (13) 

for —x 2 > R 2 . Thus, 0cm (x°,x) is independent of x° and 

0(O,x) = 0cm(-7(v-x),7[x]) = C m(O,7[x]) for |x| > R . (14) 

In the following we always consider the wave function at a; = and omit the argument for 
the relative time, for example 0(x) = 0cm (7 [ x ]) f° r ( TT4T) . We also know that the Helmholtz 
equation is satisfied : 

(V 2 + A; 2 )0 CM (x)=O for |x| > R . (15) 

We can expand 0cm (x) in terms of the spherical Bessel ji(x) and the Neumann function 

ni(x) as 



00 / 
-•cm 



00 I 

*)=Y.H M*0 • Y im{^) (ai(fc) • ji{k\x\) + A(*0 • ^(^|x|)) , (16) 



1=0 m=-l 



with some constant bi m (k), where Q is the spherical coordinate for x. The conventions of 



ji(x), ri[(x) and Yi m {Vt) agree with those of Ref. [13J. The coefficients ai(k) and fii(k) yield 
the scattering phase shift as tan^(fc) = /3i(k) / 'ai(k) . 

Next we consider a wave function on a periodic box L 3 in the Euclidian space ( on the 
lattice ), which is defined by 

vJ/ L ( Xl ,r;x 2 ,r) = (0|7r + (x l7 t)tt+(x 2 , r)|7r+7r + ; E, P) , (17) 

where |7r + 7r + ; E, P) is an eigenstate with the total energy E and momentum P on the lattice. 
The wave function is periodic with respect to the position of either of the pions : 

^ L (x 1 + n 1 L,r;x 2 + n 2 L,r) = ^ L (x x , r; x 2 , r) for m, n 2 G Z 3 . (18) 

We introduce a relative wave function like as in the Minkowski space by 

L (x) = ^( Xl ,r;x 2 ,r)-e^- iP - x , (19) 



with x = xi - x 2 and X = (xi + x 2 )/2. (JTHD and (JTSJ) yield 

L ( X ) = (_l)V(2-)P-n . ^L( x + nL ) for n e Z 3 _ 



(20) 



In the derivation of the finite size formula, it is assumed that the two-pion interaction 
range R is smaller than one-half the lattice extent, so that the boundary condition does 
not distort the shape of the two-pion interaction. With this assumption, L (x) satisfies the 
same equation as for 0(x) : 



V 2 



P ■ V) 2 /E 2 + k 2 L (x) =0 for R < | 7 [x]| < L/2 



(21) 



The finite size formula is given by solving (|21|) under the boundary condition (12171) . It is 
convenient to change the variable from x to y = 7[x] and define a new function by 



0cW(y) = L (x) , 

to solve the differential equation (12T1) . 0^ M (x) satisfies 

" 2 ]:2] aL f -'^ - for R < |x| < L/2 , 

.l) i /(2*)P-= . 0^ M (x + 7 [n]L) for n G Z 3 



v 2 + r i0cm^ x ' 



y CM> 



(22) 

(23) 
(24) 



Note that ( 1221) is not the Lorentz transformation on the lattice, but it is merely a definition 
of the function with the change of the variable from x to y = 7[x] to solve (l2Tj) . 
The general solution of (123]) under (|24|) can be written by 



oo I 



-'cm 



v x ) = /J y] v im(k) ■ Gi m (x; k) 



(25) 



1=0 m=-l 

with some constant Vi m (k). G; m (x; fc) is given from the periodic Green function : 



G(x; fc) 

r P = 



,7\3 Z_^ „2 _ 



a «qx 



7L 3 ' — ' q^ 



A; 2 



q 



■in 



q= ,~ L |— n + P/2 



nG Z 3 



(26) 

(27) 



as 



G, TO (x;fc)=;y, m (V)G?(x;AO , (28) 

where ^ m (x) = |x| z • Yi m (Q) with the spherical coordinate Q for x. The expansion of G(x; fc) 
in terms of j/(x) and n^(a;) is given by 
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00 z 



G(x; fc) = — n (/c|x|) + ^ ]T a m (fc; 1) • Y lm (^)ji(k\^\) . 



47T 

(=0 m=-J 

The function gi m (k; 1) is an analytic continuation of 



(29) 



gira(k;z) 



Air 
k ) 7L 3 



E 3Uq) " (q 2 - fc 



2\-2 



(30) 



qer F 



which is defined for Re(z) > (/ + 3)/2, where the region Tp is defined by ( |27|) . The explicit 
expansion for £q m (x; k) with general I and m is not needed. Note, however, that G/ m (x; k) 
contains j//(A;|x|) for a whole range of /' and n^(A;|x|) with only /' = / as known from (|28]) 
and (129]). 

The wave function in the Euclidian space is related to that in the Minkowski space by 
the analytic continuation of the relative time ix° — > x 4 . It seems that L (x) = 0(x) at 
x° = x 4 = and 

0cm( x ) = 0cm(x) , (31) 

from (TT4|) and (I221) . But this is not true, because the degeneracies of the energy eigenstate 
in the infinite and the finite volume are different. (|3~T1) should be changed to 



oo / 

-"CM ~" 



x ) = E E C *™(*) ■ ^Wff(M) , (32) 



Z=0 m=-Z 



with some constant Ci m (k), where f2 is the spherical coordinate for x. 0qm (l x l) * s the ^ m 
component of 0cm (x) defined by 

0CM(|x|) = |dfiF^(fi)0 C M(x). (33) 

Substituting ({TBI) and (1251) into (1321) . we obtain 



oo £ 

'CM 

Z=0 m=-Z 
oo I 



x) = ^J ^J Vimik) ■ G/ m (x; fc) 

Z=0 m=-l 
oo 2 

J] ^ D lrn (k)-Y lm (n)U(k)-ji{k\x\)+Pi(k)-ni{k\x\)) , (34) 



/=0 m=-J 

where D im (k) = Ci m (k)bi m (k) . 

In the present work we consider only the two-pion state in the Af representation of the 
rotational group on the lattice, which equals to 5-wave up to angular momentum 1 = 2. If the 
scattering phase shift for I > 2 is very small in the energy range under consideration, j3i{k) ~ 
for I > 2 in (1341 . This also means that vi m (k) ~ for I > 2, because G; m (x; k) contains 
ni(k\x\) as discussed before. This expectation is supported by our numerical simulation as 
shown later. 

We use the expansion form of G(x; k) in ( 1291) to determine the allowed values of k, Di m (k) 
and v 00 (k) in (134"|) . Comparing the coefficients of j (^l x |) and n (/c|x|) of both lines of (1341) . 
we find 

A)o(fc) • ao(fc) = v 00 (k) ■ - 7 =g 00 (k; 1) , (35) 

V47T 



D 00 (k) ■ (3 (k) = v 00 (k)-— , (36) 

47T 



where goo(k', 1) is defined from (!30|) . Finally we obtain the finite size formula by taking the 
ratio of (1551) and (1551) : 

Mk) ~ ^W) ~ — 9oo{k; 1} • (37) 

The other components of (1341) give only information for Di m (k). In the case of P = 0, this 
formula turns to be the formula presented by Luscher in Ref. [7J. 

III. METHOD OF CALCULATIONS 
A. Calculation of wave function 

In the present work we consider the ground state of the 1 = 2 5-wave two-pion system 
with the total momentum P = and P = (2ir/L)e x . When the interaction between two 
pions is turned off, the scattering momentum k defined by flHj) takes 

k = for P = , 

k = ir/(jL) for P = (2n/L)e x . (38) 

These values of k are changed by the two-pion interaction. 

In order to calculate the wave function we construct the correlation function : 

F_(x,r,r s ) = (0|fi(x,r)O(P,r s )|0) . (39) 

The operator fi(x, r) is defined by 

n ^ r ) = E7/E e4P ' X7r+ ( x+ ^ x ]' r ) 7r+ ( x ' r ) ' ( 4 °) 

R x 

where vr + (x, r) is an interpolating operator for 7r + at the position (x, r). The operation R 
represents an element of the cubic group (Oh) for P = and the tetragonal group (D^) 
for P = (27i/L)e x . The summation over R projects out Af representation of these groups, 
which equals to the 5-wave state up to the angular momentum I = 4 for P = and / = 2 
for P = (2n/L)e x . 

The operator ft(P, T s ) in (1391) is defined by 

-. N R 

Q(P,r s ) = — X)[T + (P^.;^i)T + (0,T.;i7 i )] t , (41) 



N R 



where 



7r+(P,T.;&) = i^ e ' p %r s )(;(x)] 75 [5] ! i(y 1 r s )( J (y)] . (42) 

x y 

The operator 7r + (P, t s ; rjj) is defined as 7r + (P, r s ; £j) by changing £j(x) to rjj(x). The func- 
tions ^j(x) and r]j('x) are U(l) noise whose property is 

1 N R 

J im F^^ x ^ y) = * 3(x - y) ' (43) 

3=1 



In the present work we set Nr = 2. 

Neglecting the contributions from excited states, we can extract the wave function by 

^(x) = p't'^ , (44) 

-fV7rl x Cb r 0) T s) 

in the large To region, introducing a reference position Xo- We note that L (x) at all positions 
x are not independent. The number of independent positions is (L + 2)(L + 4)(L + 6)/48 
for P = and (L + 2) 2 (L + 4)/16 for P = (2n/L)e x , owing to the boundary condition (J20J) 
and the rotational symmetry on the lattice, L (x) = <fr L (R{x]). 

In the present work we attempt to extract the energy of the two-pion system and the 
scattering phase shift from the wave function. For a comparison, we also evaluate them from 
the two-pion time correlator : 

G n7T (r, T a ) = -^22 -^( x ' t, T s ) , (45) 

x 

as done in the previous works of the scattering phase shift. We also calculate the time 
correlator for the pion and the p meson, 



GAr,r s ) = ^EjlE( l 7r+ ( x ' r )( 7r+ ( O ' T -^)) t | } , (46) 

R 3=1 x 

G P (r,r s ) = —J2j^J2 E (0|Pa(x,r)(p a (O,r s ;O)) |0> , (47) 



n j=\ ' x o=l,2,3 



which are used to extract the masses. In fT4T|) p a (x, r) is an interpolating operator for the p 
meson with the polarization a at the position (x, r) and p a (0, r s ; £,) is given by substituting 
7 a for 7 5 in 



B. Simulation parameters 

Our simulation is carried out in the quenched approximation with the plaquette gauge 
action at j3 = 5.9. Gauge configurations are generated with the 5-hit pseudo-heat-bath 
algorithm and the over-relaxation algorithm mixed in the ratio of 1 : 4. This combination is 
called a sweep and the physical quantities are measured every 200 sweeps after 2000 sweeps 
for the thermalization. 

We use the improved Wilson action for quarks [lij . The clover coefficient Csw is chosen 
to the mean-field improved value defined by 

Csw = (Un)- 3/4 = (1 - 2/(3)-"" = 1.364 , (48) 

where (Uu) is the lxl Wilson loop which is evaluated in one-loop perturbation theory. 
The quark propagators are calculated with the Dirichlet boundary condition imposed in the 



time direction and the periodic boundary condition in the spatial one. The source operator 
fi(P, T s ) in (1391) is set at r s = 20 to avoid effects from the time boundary. 

The lattice cutoff is estimated as 1/a = 1.631(45) GeV (a = 0.1208(33) fm) from m p . The 
lattice size is 32 3 x 120, which corresponds to the spatial extent 3.87 fm in the physical unit. 
We choose three quark masses to give m^ = 0.420, 0.488 and 0.587 GeV . The numbers of 
configurations are 400, 212 and 212 for each quark mass, which are generated independently. 
The masses of the pion and p meson calculated from the time correlator are listed in Table [H 

IV. RESULTS 

A. Wave functions 

We obtain the wave function L (x) by f|4*4|) and transform it to 0^ M (x) = L (7 _1 [x]). In 
the transformation we calculate the Lorentz boost factor by 7 = \/E 2 — P 2 with the energy E 
extracted from the two-pion time correlator. In Fig. [I] we show 0q M (x) at m^ = 0.420 GeV 
for the total momentum P = and P = (2n / 'L)e x . Here we set the reference position 
xo = (7, 5, 2) and tq — r s = 40 in ( 144"}) . confirming that the wave function does not depend on 
the choice of r for r — r s > 32. In the figure the left and right panels for each momentum 
show the wave functions on xy-pla.ne at z = and yz-pl&ne at x — 0. While the boundary 
conditions for P = are the periodic conditions in all directions, those for P = (2ir/L)e x 
are anti-periodic in the x direction and periodic in the other directions as known from ( !24l ). 
These features are clearly shown in the figure. 

We consider the two-pion interaction from the ratio : 

v ° M = IS? ^( v2 - (p ■ v)2/£2 > i(x) • (49) 

where y = 7[x] and V 2 , is the Laplacian with respect to y. Away from the two-pion 
interaction range, we expect that V^ M (x) is independent of x and equals to — k 2 . In the 
calculation of the ratio, we rewrite (J4"9l by 

^>^(^S + lrSy< x >- < 5o > 

and adopt the naive numerical derivative, 

' /(x) = /(x + i) + /(x-i)-2/(x), (51) 



dx: o 

where the Lorentz boost factor 7 is calculated with E extracted from the two-pion time 
correlator. In FigfJ] V^ M (x) is plotted for the same parameters as for Fig. [lj We find that 
the ratio is almost constant for |x| > 10 both for P = and P = (27r/L)e x . We also observe 
a strong repulsive interaction near the origin consistent with the negative scattering phase 
of the 1 = 2 two-pion system. 



Here, we note a physical meaning of Vqm(x) in (H9l) . It seems that Vqm(x) equals to the 
corresponding ratio in the Minkowski space defined by 

VVcm(x) 

which approximately takes the potential of the two-pion interaction in the non-relativistic 
limit. This, however, is not true for P 7^ in the interaction region |x| < R. Ignoring the 
difference of the degeneracy of the energy eigenstate in the infinite and the finite volume, 
the wave function in the Euclidian space, which is calculated in the present work, is related 
to that in the Minkowski space by the analytic continuation as </> L (x) = 0(x) at the relative 
time x° = 0. But a relation between 0^m( x ) an d 0cm (x) is 



lL ^r„n _ ml 

; CM 



(7[x]) = L (x) = 0(x) = C m(-7(v • x), 7 [x]) , (53) 



which is not a relation at a; = 0. Thus, there is no simple relation between V^ M (x) in (149p 
and Vcm( x ) i n f|52|) for P ^0. Only for the region |x| > R, we give a simple relation : 
Vqm(x) = Vcm(x) = —A; 2 , due to the x°-independence of the wave function, 0cm(^°,x) = 
0cm(O,x). 

We now consider the two-pion interaction range R. In the quantum field theory the wave 
function does not strictly satisfy the Helmholtz equation (1231) . even for the large |x| region. 
Hence, with k obtained from the two-pion time correlator, V^ M (x) + k 2 shows a small tail at 
large |x|. We may take the wave function as satisfying the Helmholtz equation, HV ( j M ('x) + k 2 
is sufficiently small compared with k 2 . In the present work we take an operational definition 
of the range R as the scale, where 

EWx) = p = , 2#m(x) • («) 

is small enough compared with the statistical error. With this definition we expect that the 
systematic error of the scattering phase shift from the interaction tail is smaller than the 
statistical error of the scattering phase shift. 

We show Uq M (x.) as the function of |x| in Fig. [3] We find R < L/2(= 16) in all cases. 
This means that the necessary condition for the finite size formula is satisfied on the 32 3 
lattice for our range of the quark mass m n = 0.420 — 0.587 GeV and the momentum k ~ 
— 7r/(7L)(~ 0.16 GeV) with the current statistics of the simulations. 

B. Scattering phase shift from wave function 

We attempt to calculate the scattering phase shift by substituting the momentum k ob- 
tained with three methods into the finite size formula ( 1371) . First method is the conventional 
one, where we calculate the momentum k by ([6]) with the energy of the two-pion system E 
extracted from the two-pion time correlator. Our results of k and the scattering phase shift 
#o(A;) are tabulated in Table HTl for P = and Table HTT1 for P = (2n/L)e x (labeled "from 
T"). 

10 



In the second method, we extract k by fitting the wave function L (x) to the Green 
function (1261) with the form, 

x)=C-G(7[x];A0, (55) 



L.L, 



taking k and an overall constant C as free parameters. This method was introduced by the 
CP-PACS collaboration in Ref. [9( for the evaluation of the scattering length, where they 
found that the statistical error can be reduced with this method. The numerical evaluation 
of G(7[x]; k) and an explicit procedure of the fitting are discussed in Appendix|A] We choose 
the fitting range |x| > x m with x m tabulated in Table HT1 and IIHI An example of the fitting 
is shown in Fig. H] at m n = 0.420 GeV, where the data points are shown with open circles 
and the values for the fits with cross symbols. We find that the fit works well both for P = 
and P = (2ir/L)e x . These mean that the contribution of G/ m (x; k) with / > 2 are negligible 
as expected. Our results of k and So(k) are tabulated in Table HT1 and IIHI (labeled "from 0"). 
These are consistent with those of the conventional method ("from T"). We also find that 
the statistical errors is reduced. Thus the method with the wave function is also efficient for 
the evaluation of the scattering phase shift. 

As discussed in Sec. \IV A\ the ratio Vq M ('x) in ( 1491) equals to — k 2 away from the interaction 
range. For simplicity we write the ratio by 

K(x; 7) . l&ttH) = ^ (JL*, + J + Q ^( X ) , (56) 

with (l50l) . We extract k by minimizing the chi-square : 

|x|>a; m 

taking k as a free parameter, where SV(x.; 7) is the statistical error at fixed 7. This method 
is an extension of that introduced by the CP-PACS collaboration in Ref. J9(] for P = to 
the system with P/0. As known in fl56l) . V(x; 7) depends on k through 7. This fact makes 
the analysis difficult. We use a similar procedure as for fitting the wave function (see in 
Appendix [B]). We choose the same fitting range x m for fitting the wave function. Our results 
of k and So(k) are tabulated in Table HI1 and IIHI (labeled "from V"). As shown in the tables, 
the results are consistent with those from the two-pion time correlator ( "from T" ) and fitting 
the wave function ("from <p n ). We also find that the statistical error is significantly reduced. 

C. Scattering length and phase shift at physical quark mass 

For P = 0, the momentum k is very small as shown in Table [Til Thus the scattering 

length a can be calculated by 

tan 5 (k) 
a = . (58) 

In Fig. Owe plot the scattering length aQ/m n obtained from three methods discussed in the 
previous section. They are consistent, but the statistical errors are significantly reduced by 

11 



using the two methods with the wave function ("from 0" and "from V"). We carry out the 
chiral extrapolation with the fit form, a^jm^ = A + B ■ m 2 and obtain 

a / m ^ = -2.211(77) GeV" 2 , (59) 

in the chiral limit, for the data of fitting V(x, 7) ("from V"). The result of chiral ex- 
trapolation is also plotted in the Fig. [5j We refer that the prediction from CHPT is 



a /m w = -2.265(51) GeV" 2 15 



We obtain the scattering phase shift at the physical quark mass for the various momenta 
k from the scattering amplitude defined by 

tan<5o(&) 



Aim*, k) = iau ^ w ^ + ^ , (6°) 

which is normalized as Aim^, k) = a^vn^ at k — 0. We show the results in Table HT1 and IHH 
and plot in Fig. [61 We fit these amplitudes with a fitting form : 

A(m n , k) = A w ■ ml + A 20 ■ m\ + A 01 ■ k 2 + A n ■ mlk 2 . (61) 

The results of the fitting are also plotted in Fig. [6] and summarized in Table IIVI We find 
that Aio(= do/ ni-n) is consistent with (1591) . 

Finally, we estimate the scattering phase shift $o(k) at the physical quark mass by putting 
m-x = 0.140 GeV in flBTj). In Fig. [7] our results are plotted as a function of k 2 and compared 
with the experiments [l6j. Our results are larger than the experiments. We consider that 
this disagreement arises from the effect of the finite lattice spacing. We must leave this issue 
to studies in the future. 



V. CONCLUSION 

We have calculated the two-pion wave function for the ground state of the 1 = 2 S-w&ve 
two-pion system both for P = and P = {2ir/L)e x . We have investigated the validity of 
the necessary condition for the finite size formula and found that it is satisfied on the 32 3 
lattice for the quark mass range m n = 0.420 — 0.587 GeV and the scattering momentum 
k < 0.16 GeV. We have also found that the scattering phase shift can be extracted from the 
wave function with a smaller statistical error than from the two-pion time correlator, which 
have been used in the studies to date. 

An implication of the present work is the feasibility to calculate the decay width of 
the p meson through studies of the 1 = 1 two-pion system. While the evaluation of the 
disconnected diagrams with a good precision has been a computational problem, our method, 
investigating the scattering system from the wave function in which the energy eigenvalue 
is extracted from the wave function at a single time slice, could lend a tactics that can be 
used to evaluate such complicated diagrams with a modest cost. 
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APPENDIX A: FITTING WAVE FUNCTION TO GREEN FUNCTION 

We introduce 

P=(r l [ p /2]) 2 ) (A1) 

and expand the Green function around k 2 = k 2 as 

oo 

G( 7 [x]; k) = J> 2 - Fy- 1 • F( 7 [x];j, k) . (A2) 

i=o 



(A3) 



The coefficient F(x; j, fc) is defined by 




F(x;0,fc) = 


1 V e iqx 
7 L 3 ^ 

' q 2 =fc 2 


and 


X- ! 



for j > 1, where q e Tp (defined by (p7|) ). We evaluate F(x; j, fc) for j > 1 from 

( 7 L 3 )-F(x;j,fc) 



(A4) 



e jqx 



+ TrT)\ (t) ~ 3 jf d " " 2J " e< """ ' /(xi "> • ( A5) 



where 



/(x; 77) = X] (-l) L /^ p - n ■ e -(( x +^ n ]i)/(2r,)) (A6) 

n€Z 3 

5]) is obtained by the same technique discussed in Ref. [9|. F( 7 [x];j, fc) in the (1A2|) is 



depend on k through 7 . This fact makes the fitting of the wave function difficult. In the 
present work, we adopt the following procedure: 



[1] We calculate 7 o = Ej\/E 2 — P 2 with the energy E extracted from the two-pion time 
correlator and take it as the initial value of 7 . 

[2] We fix k and F( 7 [x]; j, k) in the (1A2I) with given 7 and fit the wave function, taking 
k as the free parameter. 



[3] We update 7 by 7 = y/1 + P 2 /[A(ml + k 2 )] with k obtained in (2J. 

[4] We iterate the procedure [2] - [3] , until the value of k becomes to be stable. 

It is expected that k rapidly converges in this procedure, because the dependence of 7 on k 
is very small. We confirm that k is stable within the single precision after five iterations for 
all our simulation parameters. We carry out above procedure for the each jackknife bin and 
estimate the statistical error of k by the jackknife method. 
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APPENDIX B: FITTING V(x; 7) 

The fit of V(x;7) is performed by the similar procedure as for fitting the wave function 
discussed in Appendix [A] : 



[1] We calculate 70 = Ej\/E 2 — P 2 with the energy E extracted from the two-pion time 
correlator and take it as the initial value of 7. 

[2] We fix V^(x; 7) with given 7 and search k that gives the least number of x 2 (/c) in (I5TJ) . 



[3] We update 7 by 7 = a/1 + P 2 /[4(m2 + k 2 )} with k obtained in (2J. 

[4] We iterate the procedure [2] - [3] , until the value of k becomes to be stable. 

It is expected that k in this procedure rapidly converges, because the dependence of 7 on k 
is very small. We confirm that k is stable within the single precision after five iterations as 
for the fitting the wave function. We carry out above procedure for the each jackknife bin 
and estimate the statistical error of k by the jackknife method. 
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FIGURES 



(a) P = 





(b) P = (2ir/L)e x 





FIG. 1: Two-pion wave functions 0q M (x) at i% = 0.420 GeV for the total momentum P = and 
P = (2ir/L)e x . The left and right panels for each momentum show 0^ M (x) on xy-plane at z = 
and yz-plane at x = 0. 
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(a) P = 





(b) P = (27r/L)e a 





FIG. 2: V^ M (x) at m^ = 0.420 GeV for the total momentum P = and P = (2w/L)e x . The left 
and right panels for each momentum show Vq M (x.) on xy-plane at z = and yz-plane at x = 0. 
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(a) P = 



(b) P = (27r/L)e a 
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FIG. 3: Uq M (x.) as the function of |x| for the total momentum P = and P = (2n/L)e x . 
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(a) P = 
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Ixl 



(b) P = (2vr/L)e a 
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8.0 10.0 12.0 14.0 16.0 18.0 20.0 22.0 24.0 26.0 28.0 

|x| 

FIG. 4: Results of fitting the wave function L (x) to the Green function at m^ = 0.420 GeV for 
the total momentum P = and P = (2ir/L)e x . Open circles refer to the data points and cross 
symbols to the results of the fitting. 
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0.2 0.3 

m/ (GeV 2 ) 



0.5 



FIG. 5: Scattering length ao/m n obtained from the two-pion time correlator ("from T"), fitting 
the wave function ("from (/>") and fitting V(x;j) ("from V"). The result of chiral extrapolation 
for the data of "from V" is also plotted. 
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FIG. 6: Scattering amplitude ^4(m T , k) fore three quark masses obtained from the two-pion time 
correlator ("from T"), fitting the wave function ("from </>") and fitting y(x;7) ("from V"). The 
fit curves for the data of "from V" are also plotted. 
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FIG. 7: Results of the scattering phase shift at the physical pion mass and the experiments in 
Ref. EF 
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TABLES 



0.1413 
0.1410 
0.1405 



jVconf. 
400 

212 
212 



m* (GeV) 
0.41964(66) 
0.48829(95) 
0.58749(74) 



m p (GeV) 
0.837(12) 
0.8598(95) 
0.9007(67) 



m n /m p 
0.5012(71) 
0.5679(62) 
0.6522(48) 



TABLE I: Masses of the pion and the p meson calculated from the time correlator. 



K 




0.1413 


0.1410 


0.1405 


m n (GeV) 




0.41964(66) 


0.48829(95) 


0.58749(74) 


%m 




13.0 


14.0 


12.0 


k (GeV) 












from T 


0.0383(14) 


0.0360(19) 


0.0396(14) 




from (f> 


0.03760(91) 


0.0363(12) 


0.03810(92) 




from V 


0.03762(46) 


0.03750(57) 


0.03875(33) 


6 (k) (deg.) 












from T 


-1.71(18) 


-1.44(22) 


-1.88(18) 




from (f> 


-1.63(11) 


-1.48(14) 


-1.70(11) 




from V 


-1.636(56) 


-1.621(69) 


-1.776(42) 


ao/m-K (GeV -2 ) 












from T 


-1.86(12) 


-1.43(14) 


-1.415(87) 




from (f> 


-1.808(79) 


-1.460(88) 


-1.322(58) 




from V 


-1.809(40) 


-1.546(42) 


-1.363(21) 


A(m n , k) 












from T 


-0.330(22) 


-0.343(34) 


-0.489(30) 




from (j) 


-0.320(14) 


-0.349(21) 


-0.457(20) 




from y 


-0.3198(70) 


-0.370(10) 


-0.4709(72) 



TABLE II: Results for the total momentum P = from two-pion time correlator ("from T"] 
fitting wave function ("from (j)") and fitting V(x;7) ("from V"). 
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K 




0.1413 


0.1410 


0.140 


m* (GeV) 




0.41964(66) 


0.48829(95) 


0.58749(74) 


%m 




10.0 


11.0 


12.0 


k (GeV) 












from T 


0.1610(14) 


0.1594(21) 


0.1634(12) 




from (/> 


0.16015(96) 


0.16117(82) 


0.16369(54) 




from V 


0.15956(40) 


0.16197(38) 


0.16394(28) 


7 


from T 


1.06159(22) 


1.04749(23) 


1.033921(92) 




from 4> 


1.06168(19) 


1.04739(17) 


1.033944(82) 




from y 


1.06173(17) 


1.04734(16) 


1.033937(82) 


So(k) (deg.) 












from T 


-8.4(12) 


-5.5(17) 


-7.20(98) 




from <j> 


-7.69(78) 


-6.92(69) 


-7.46(46) 




from V 


-7.21(33) 


-7.59(31) 


-7.67(24) 


A(m n , k) 












from T 


-0.412(55) 


-0.308(94) 


-0.472(62) 




from </> 


-0.379(37) 


-0.387(37) 


-0.488(29) 




from y 


-0.356(15) 


-0.423(17) 


-0.501(15) 



TABLE III: Results for the total momentum P = (2Tr/L)e x from two-pion time correlator ("from 
T"), fitting wave function ("from 0") and fitting V(x;7) ("from V"). 



Aio (GeV" 



A 20 (GeV" 



An (GeV- 2 ) 



A u (GeV~ 



from T 
from (j) 
from V 



-2.10(27) 
-2.17(17) 
-2.187(80) 



2.02(92) 
2.52(59) 
2.42(26) 



-7.7(56) 
-3.2(34) 
-2.3(15) 



26(22) 

6(12) 

3.2(56) 



TABLE IV: Results of the coefficients Aij in ([61 
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